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Abstract* Species diversity and abundance of fruit-feeding nymphalid 
butterflies were studied in an Amazon rain forest fragment. Butterflies 
were caught in baited traps in twelve areas, selected to sample a gradient 
of increasing disturbance. Measurements of six parameters of vegetation 
structure were also taken to estimate the disturbance. A total of 90 but¬ 
terfly species were trapped. The greatest alpha diversities were found at 
the edge of the forest and in areas of intermediate disturbance. Canoni¬ 
cal Correlation Analysis (CCA) showed that the composition of the spe¬ 
cies assemblages of nymphalids was related to vegetation structure vari¬ 
ables, especially girth at breast height and number of tree morpho-spe- 
cies. The butterfly fauna appeared more similar in forested areas than in 
the disturbed ones. Some species were suggested as habitat indicators and 
the value of this guild of fruit-feeding butterflies in conservation programs 
is discussed. 

Key Worbs: Butterflies, nymphalids, diversity, community structure, distur¬ 
bance, rain forest fragment, direct gradient analysis, Brazil. 

Introduction 

One of the main objectives of community ecology is the synthesis of the 
roles of physical and biological factors that determine species abundance 
and distribution within and among natural communities. After MacArthur 
and MacArthur (1961) found a relationship between bird diversity and high 
vegetation diversity, ecologists verified that habitat complexity is an impor¬ 
tant factor for the structuring of local communities. Habitats that are struc¬ 
turally more complex and heterogeneous offer more niches, and therefore 
support a greater number of species (spatial heterogeneity, Fianka 1966). 
In addition this idea, Connell (1978) suggested that high diversity in tropi¬ 
cal forests is maintained by disturbances, such as tree falls. Considering such 
dynamics, the forest can be seen as a mosaic of gaps in different succes- 
sional stages, with different local communities, and a high regional diver¬ 
sity. 

Although biogeographic and historical conditions are extremely impor¬ 
tant factors in structuring communities (Slansky 1972, Ricklefs 1987, Leps 
and Spitzer 1990, Brown 1982, Brown 1991, Thomas 1991, Gaston 1996), 
local factors also affect local butterfly diversity (Emmel and Leek 1969, 
Montesinos 1985, DeVries 1994, Kitahara and Fuji! 1994, Sparks and Par- 

Paper submitted 7 July 1997; revised manuscript accepted 25 February 1999. 


30 


j. Res . Lepid. 


ish 1995). Many studies have shown that tropical butterfly communities also 
respond to physical factors of the habitat, such as topography, stratification, 
gaps, edges, urbanization and habitat disturbances (Ruszczyk 1986, DeVries 
1988, Raguso and Llorente-Bousquets 1990, Brown 1991, Hill et al. 1992, 
Pinheiro and Ortiz, 1992, Hill et al. 1995, Spitzer et al. 1997). As such, 
multivariate analysis has proven to be an important tool when investigating 
the relationships between species assemblages and environmental variables 
(Leps and Sptizer 1990, Kremen 1992, Ramos 1992, Vaisanen 1992, Spitzer 
et al. 1993, Spitzer et al. 1997, Blair and Launer 1997). 

The objectives of this study are (1) to measure neotropical nymphalid 
butterfly diversity along a gradient of disturbance, (2) to explore the spe¬ 
cies-environment relationships through a direct gradient analysis and (3) 
identify the most important butterfly species and vegetation variables, which 
could be used in conservation monitoring programs. 

Methods 

Study site. This study was conducted in a forest fragment at the boundary of the 
eastern Amazon (5°0DS, 47°32AV; 260 m), a region where the natural landscape 
has been greatly modified by human activity. The study site was about 50 km north 
of the transition to Cerrado. The fragment has about 1,000 ha of primary forest 
with several levels of disturbance, surrounded by secondary forest in several 
succesional stages, eucalyptus monocultures and cattle pasture. I selected 12 sample 
units (SUs) throughout a disturbance gradient: forest understory (FIJI, FU2, FU3), 
forest roads (FR1, FR2, FR3), edge (EDO), highly disturbed forest understory 
(DFU), highly disturbed forest road (DFR), 4-year-old secondary forest (SF4), 2- 
year-old secondary forests (SF2) and Eucalyptus pellita monoculture (EUC). 

Data collection. I made lepidoptera collections between June 1990 andjuly 1991. 
For each of 12 SUs, three fruit-baited traps were set in line, suspended 1.0-1.7 meters 
above the ground, and 25 meters apart from one another. For each collection, the 
traps were visited for 14 consecutive days. The banana and sugar cane bait was kept 
moist for the duration of the trapping period. The disturbance level of each SU 
was estimated using vegetation parameters obtained through the point-centered 
method (Miiller-Dumbois and Ellemberg 1974), with 21 quartered points estab¬ 
lished per SU, only for trees up to 20 cm of circumference at breast high. This 
method was chosen for its simplicity and common use in phytosociological surveys. 
The following vegetation variables were used: average girth at breast hight (GBH); 
estimated average tree height (THG); number of tree morpho-species (NMS) estimated 
by rind and leaf characteristics, with the help of a local guide; tree density within 100 
m 2 (DEN); average horizontal cover (HOC), estimated at each sample point by an 
observation made on a 50 cm square carton held 10 m from the observer in each 
quarter. Cover was estimated to be within one of four categories (0-25%, 25-50%, 
50-75% and 75-100% vegetation cover); average vertical cover (VEC), estimated by 
the four previously mentioned vegetation cover categories, apllying a 10 cm square 
frame held at a distance of 60 cm from the observer at an angle of approximately 
20° in relation to zenith. 
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Table 1. Alpha diversity of fruit-feeding butterfly species in twelve sample units 
of an Amazonian forest fragment. 


Forest Forest Edge Disturbed Secondary Eucalvp 

understory road forest forest monoc 



FU1 

FU2 

FU3 

FRl 

FR2 

FR3 

EDG DFU 

DFR 

SF4 

SF2 

EUC 

Total 

Number of individuals (N) 63 

43 

97 

114 

111 

106 

334 

267 

571 

490 

744 

604 

3544 

Nymphalinae richness 

3 

5 

6 

11 

13 

12 

15 

16 

19 

17 

16 

12 

29 

Satyrinae richness 

13 

11 

15 

14 

16 

13 

27 

14 

20 

19 

20 

17 

41 

Brassolinae richness 

1 

1 

9 

0 

1 

1 

5 

9 

4 

3 

2 

9 

6 

Charaxinae richness 

0 

2 

3 

2 

3 

6 

6 

6 

6 

7 

3 

3 

13 

Total species richness (S) 

18 

20 

27 

28 

34 

33 

54 

39 

50 

47 

42 

34 

90 


Species diversity (H') 2.27 2.76 2.60 2.79 3.08 3.06 3.38 3.06 3.19 3.00 2.28 2.09 3.26 


Data analysis. The butterfly alpha diversity of each SU was quantified by the spe¬ 
cies richness (S) and Shannon-Wiener index (H’). To evaluate environmental ef¬ 
fects on the butterfly community I ran a Canonical Correspondence Analysis (CCA), 
using the program CANOCO (TerBraak 1988). The vegetation parameter estimates 
of habitat disturbance were used as environmental variables in the CCA. The vari¬ 
able tree height was removed from the analysis due to its high value of inflation, 
and high colinearity with the other variables. The significance of species-environ¬ 
ment relationships was tested using a Monte Carlo test. 

Results 

A total sample effort of 2,016 trap days (=3 traps x 12 SUs x 56 days) 
resulted in 3,544 individuals collected, representing 90 species of five 
subfamilies of Nymphalidae. The five most abundant species were 
Paryphthimoides phronius, Yphthimoides spl, Yphthimoides disaffecta, 
Hermeuptychia hermes and Cissia penelope, all belonging to the subfamily 
Satyrinae, representing 45.3% of the total number of individuals col¬ 
lected. A complete list of species abundances in each SU can be seen 
in Appendix 1. The total butterfly diversity in the rain forest frag¬ 
ment sampled was FT = 3.258. 

An analysis of alpha diversity showed that edges and areas of intermedi¬ 
ate disturbance presented higher species richness and diversity (Table 1). 
Although more disturbed areas, such as eucalyptus monoculture and 2-year- 
old secondary forest had higher species richnesses than forest, they had 
lower species diversities, due to the high dominance of the Satyrinae spe¬ 
cies. This pattern was not found when other groups were considered sepa¬ 
rately: Satyrinae and Brassolinae had higher species richnesses in the dis¬ 
turbed areas, with a peak of the edge; Charaxinae and Nymphalinae had 
higher species richnesses in the road, edge and disturbed forest areas. 

The CCA ordination diagram shows the relationships between butterfly 
species, sample units and environmental variables (Figure 1). By compar¬ 
ing the arrow lengths, one may evaluate the significance of the constrain¬ 
ing vegetation variables. The arrow points roughly in the direction of the 
maximum variation in the value of the corresponding variable. The spe- 
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Figure 1. CCA ordination of the fruit-feeding butterfly communities in an Amazo¬ 
nian forest fragment, with respect to five vegetation variables (arrows). 
The species are abreviated according to the first three letters of their 
generic and the first three letters of their specific names (see Appendix 
1 for full names). The vegetation variables are: average girth at breast 
height (GBH), estimated average tree height (THG), number of tree 
morpho-species (NMS), tree density within 100 m 2 (DEN), estimated 
average horizontal cover (HOC) and estimated average vertical cover 
(VEC). The sample units are forest understory (FU1 ; FU2, FLJ3), forest 
roads (FR1, FR2, FR3), edge (EDG), highly disturbed forest understory 
(DFU), highly disturbed forest road (DFR), 4-year-old secondary forest 
(SF4), 2-year-old secondary forests (SF2) and Eucalyptus pellita monoc¬ 
ulture (EUC). 


cies far from the origin are rare and less important to analysis. The Monte 
Carlo permutation test showed that those species are significantly related 
to supplied vegetation variables (99 permutations, P < 0.01). 

In the first CCA axis, girth at breast height presented the highest abso¬ 
lute value, followed by the number of tree morpho-species, both contribut¬ 
ing to the species data fit. This axis clearly shows a gradient of disturbance 
— preserved areas with a greater richness of large bole trees on the posi¬ 
tive side, and disturbed or early succesional ones on the negative side. The 
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ordination also shows how species respond to vegetation variables: with in¬ 
vader species typical of open areas (the small Satyrinae Hermeuptychia 
hermes , Yphthimoides spl , Y. disaffecta , Cissia penelope, Erichtodes numeria and 
Pharneuptychia pharnaces , and the Nymphalinae Hamadryas feronia., H. 
fehrua and Biblis hyperia) showing negative scores. With positive scores, 
near the origin, are the heliophyllous species of the disturbed forests, gaps, 
edges and canopy (. Eryphanis polyxena ., Hamadryas iphthime , H velutina , Mem¬ 
phis morons, Narope cyllaharus , Mca flavilla , Pareuptychia ocirrhoe , Temenis laothoe 
and Taygetes laches). The species on the right side of the diagram are typical 
of the forest understory ( Colohura dirce , Morpho achiles , Nessaea ohrinus , Taygetis 
celia, T, echo and T. virgilia) . The second axis was primarily related to tree 
density and horizontal cover, but did not form a clear gradient. The analy¬ 
sis gave a large weight to some Satyrinae species, such as Cissia penelope and 
Paryphthimoides phronius , which had large populations in high tree density 
eucalyptus monoculture. 

Discussion 

The forest edge and intermediate disturbance forest presented higher 
values of butterfly species richness and diversity. These environments, where 
intense regeneration occurs, have high productivity and maintain high 
population levels. On the other hand, the disturbance rate is high, thus 
reducing the effect of competitive exclusion. A number of ecologically based 
hypotheses have been proposed to explain patterns of species richness and 
diversity, but not all of them are mutually exclusive (Meffe and Carroll 1997). 
The productivity-disturbance hypothesis (Huston 1994) combines elements 
of several other hypotheses, proposing that the high productivity and the 
disturbance rate conditions of forest edges and gaps result in high species 
richness. The results of this study tend to agree with this hypothesis. Addi¬ 
tionally, the mixture of forest understory umbrophyllous species with open 
area heliophyllous ones raises the local diversity. The fact that edges, gaps, 
physiognomic transitions, and disturbed and secondary forests have high 
diversity has been documented many times elsewhere (Leps and Spitzer 
1990, Raguso and Llorente-Rousquets 1990, Hill et al. 1992, Pinheiro and 
Ortiz 1992, Vainsanen 1992, Spitzer et al. 1993, Spitzer et al. 1997). For 
sunloving species, drastic changes in light intensity can act as a habitat bar¬ 
rier, while edges, gaps and canopies may be treated as a continuum of sunny, 
open area. Other species, however, remain restricted to shady environments 
(DeVries 1988). 

The number of tree morpho-species was not a good predictor of the num¬ 
ber of butterfly species. Because the group is herbivorous, butterfly - host 
plant relationships have been explored (Gilbert and Smiley 1978, Erhardt 
and Thomas 1991), although not always being meaningful (Sharp et al. 1974, 
Courtney and Chew 1987, Singer and Ehrlich 1991). Besides that, the sam¬ 
pling considered only trees, and did not consider bushes, herbs and lianas 
that are host plants of several butterflies (DeVries 1987). Although not re¬ 
lated with butterfly alpha diversity, the number of tree rnorpho-species was 


34 


J. Res. Lepid. 


important in the formation of an environmental gradient, and is related to 
habitat disturbance. 

The ordination diagram shows that the forest assemblages are homoge¬ 
neous, but among disturbed areas there are great variations in species com¬ 
position and abundance. A larger constancy of forest communities has been 
verified for neotropical satyrins (Brown 1991), and butterflies of South-east¬ 
ern Asia (Leps and Spitzer 1990), but the opposite pattern was found for 
Notodontidae and Arctiinae in the same forest fragment (Dubois 1993). 

Although other butterfly groups such as Ithominae and Heliconinae have 
been suggested as more efficient indicators (Brown 1991, Beccaloni and 
Gaston 1995), the use of the fruit-feeding Nymphalidae has its utility in 
conservation programs (Daily and Ehrlich 1995). The results of the CCA 
ordination show that butterfly communities have a significant relationship 
with vegetation variables, and suggest the use of this assemblage as an ap¬ 
propriate indicator of habitat heterogeneity over this spatial scale. CCA can 
be used to match a species assemblage to environmental factors for which 
it is a good indicator, and select a subset of species as indicators for more 
intensive monitoring (Kremen 1992). Since rare species have little weight 
in the analysis (TerBraak 1988), common species, and not rare ones, should 
be selected from this guild to be used as indicators. Thus, Hermeuptychia 
hermes , Yphthimoides spl, Y. disaffecta, Cissia penelope , Erichthodes numeria , 
Phameuptychia pharnaces, Hamadryas feronia, H. februa and Biblis hyperia may 
serve as indicators of disturbed environments, and Colobura dirce , Morpho 
achiles, Nessaea obrinus , Taygetis celia, T. echo and T. virgilia indicators of more 
preserved environments. 

The use of higher taxa for biodiversity measurements (Williams and 
Gaston 1994) can be an important management tool for situations where 
taxonomic identification at the species level is difficult. For the same data 
set, counting only subfamily abundance, Ramos (1992) obtained similar 
ordination patterns as when counting species abundance. Another advan¬ 
tage of this fruit-feeding guild is that it can easily be sampled with traps, 
simultaneously in several points. Using appropriate criteria and guidelines, 
as suggested by Sparrow et al. (1994), this nymphalid fauna may be an in¬ 
formative species subset for monitoring programs. 

The collection of vegetation variables was designed to be as simple as 
possible. Of course, other local habitat variables that are important for adult 
butterflies which could have been measured were not quantified. Among 
the physical and structural variables are the size of the area, topography, 
temperature, humidity, light, gaps, roosts and dormitories, and ground 
pattern. Important biological factors for adults include food and ovoposition 
site availability, predators and mimics. 

Local diversity is determined not only by local factors, but also by regional 
and historical factors (Ricklefs 1987). Aside from the limits of the local 
habitat structure, the local butterfly assemblage depends on the regional 
species pool and historical processes such as climatic changes, isolation, 
extinction and speciation. The rapid fragmentation of the Amazon rain 
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forest may be contributing to butterfly extinctions, especially larger species 
•with scarce resources - Morphinae, Brassolinae and Charaxinae (Brown 
1991). Alternatively, the vegetation structure of disturbed forest is suitable 
for sun-lovers, secondary and opportunistic species that may spread through¬ 
out the region. Some of these butterflies are common in open biomes such 
as the Cerrado. For example, Hamadryasfebrua , H . feronia , Erichthodes numeria 
and Hermeuptychia hermes are as abundant in cerrado strictu sensu of central 
Brazil (Pinheiro and Ortiz 1992) as in the disturbed areas of the fragment 
studied. 
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CHARAXINAE 

Mem ara Memphis arachne arachne (Cramer, 1777) 
Mem mor Memphis morvus morvus (Fabricius, 1775) 
Mem oen Memphis oenomais (Boisduval, 1870) 

Mem phi Memphis philumena corita (Fruhstorfer, 1916) 


Mem ryp Memphis ryphea (Geyer, 1834) 

Pre doo Prepona demophoon Hi’ibner, 1806 
Pre don Prepona demophon (Linneaus, 1758) 
Pre lae Prepona laertes (Hubner, 1811) 
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